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Abstract In this study, we have obtained the exact so- 
lutions of the Schrodinger equation for a multi-layered 
quantum antidot (MLQAD) within the effective mass 
approximation and dielectric continuum model for the 
spherical symmetry. The MLQAD is nano-structured 
semiconductor system that consists of a spherical core 
(e.g. Ga\- x Al x As) and a coated spherical shell (e.g. 
Gai-yAlyAs) as the whole anti-dot is embedded in- 
side a bulk material (e.g. GaAs). The dependence of 
the electron energy spectrum and its radial probabil- 
ity density on nano-system radius are studied. The nu- 
meric calculations and analysis of oscillator strength of 
intersubband quantum transition from the ground state 
into two first allowed excited states at the varying ra- 
dius, for both the finite and infinite confining potential 
(CP) as well as constant shell thickness, are performed. 
It is shown that, in particular, the binding energy and 
the oscillator strength of the hydrogenic impurity of a 
MLQAD behave differently from that of a single-layered 
quantum antidot (SLQAD). For a MLQAD with finite 
core and shell CPs, the state energies and the oscilla- 
tor strengths of the impurity are found to be dependent 
on the shell thickness. At the large core radius and very 
small shell thickness, our results are closer to respective 
values for a SLQAD that previously reported. 
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dot • Hydrogenic impurity • Intersubband electronic 
transitions 



Y. Naimi ■ A. R. Jafari 

Department of Physics, Lamerd Branch, Islamic Azad Uni- 
versity, Lamerd, Iran. 
Tel.: +98-782-5221100 
Fax: +98-782-5221232 
E-mail: y.naimi@iaulamerd.ac.ir 



1 Introduction 

In the recent years, low-dimensional semiconductors have 
attracted much attention, due to a wide range of their 
applications in electronic and optoelectronic devices. 
The nano-structured semiconductor systems such as quan- 
tum dots (QDs), quantum antidots (QADs), wells, well 
wires and antiwells have been central point of both the 
theoretical and the experimental researches. One of the 
most studied classes of heterostructures is the class of 
QDs. In QDs, energy levels in valance and conduction 
bands have discrete distribution and the band gap is 
larger than compared to bulk state because of quan- 
tum confinement effects. Quantum confinement effects 
occur when the size of the structure is comparable to 
the Bohr radius. In addition to QDs, in the past few 
years, a great attention has been devoted to the physics 
of the semiconductor QADs [TJ[2] ■ 

A QAD is a portion of the semiconductor material 
that its charge carriers are confined in all three spatial 
dimensions. So, QADs have electronic characteristics 
intermediate between those of bulk semiconductors and 
those of discrete molecules. The electronic properties of 
QADs are closely related to their size and shape. The 
potential application of QADs in the development of 
semiconductor optoelectronic devices has led to much 
study of their electronic characteristics [3HHIS]. 

One of the most studied problems during the last 
decade is the investigation of properties of semicon- 
ductor when it is doped with shallow donor impuri- 
ties [6l[7] . The thermal, optical and electrical properties 
of a doped semiconductor are dependent on impurity 
states. Consequently, the calculation of impurity states 
in semiconductors has received considerable attention 
[S1IS]. The first systematic study of hydrogenic impu- 
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rity states in a GaAs quantum well has been done by 
Bastard QUIT]. 

Many studies of the single-layered QD with a central 
donor impurity for different confining potential values, 
have been reported in 12 , 13 , 14 .15,16. Varshni in [17] 
calculated the energies for Is, 2p and 3d states of a 
QD by using of the variational method. The Varshni's 
results are well consistent with exact solutions that have 
reported in QIim[lim[12 ■ 

For multi-layered quantum dot (MLQD) case, con- 
sists of a spherical core (e.g. GaAs) and a coated spher- 
ical shell (e.g. Ga\- x Al x As) as the whole dot is embed- 
ded inside a bulk material (e.g. Ga\- y Al y As), Cheng- 
Ying and Dcr Sun Chuu 23 have solved the Schrbdinger 
equation and have calculated the exact ground sate en- 
ergy of a hydrogenic impurity located at the center of 
MLQD in the framework of the effective mass approx- 
imation. Their studies shown that as the dot radius 
approaches infinity, the state energies of an impurity lo- 
cated at the center of a multi-layered or a single-layered 
QD approach —l/n 2 Ry, where n is the principle quan- 
tum number. 

In this paper, we investigate the properties of the 
multi-layered quantum antidots (MLQAD). The MLQAD 
consist of a spherical core (e.g. Ga\^ x Al x As) and a 
coated spherical shell (e.g. Gai- y Al y As) as the whole 
antidot is embedded inside a bulk material (e.g. GaAs) 
with a hydrogenic impurity in the center. We obtain 
the energies and probability density functions for both 
the ground and first excited state of the MLQAD with 
different core and shell CPs. Then, by using of ener- 
gies and wave functions, the oscillator strength of in- 
tersubband transition involving the ground state to the 
first allowed excited state are computed. Our results 
show that for small core radii the oscillator strengths 
are dependent to shell thickness and core CP values but 
when the core and shell radius are several times of Bohr 
radius, the oscillator strengths are independent of the 
shell thickness and core CP. 

For many semiconductor quantum heterostructures, 
such as GaAl/Gai- x Al x As, the polarization and image 
charge effects can be significant in the multi-layered sys- 
tem if there is a large dielectric discontinuity between 
the dot and the surrounding medium. However, this is 
not the case for the GaAs/Gal x Al x As quantum sys- 
tem [21] , therefore these effects can be ignored safely in 
our calculation. Furthermore, for the sake of generality, 
the difference between the electronic effective mass in 
the antidot core, the shell material and the bulk ma- 
terial has been ignored (i.e., fi% ~ pii ~ fi^ = [i). The 
effective atomic units are used throughout the paper so 
that all energies are measured in units of the effective 
Rydberg, Ry = ^e 4 /2e 2 h 2 , and distances are in units 



of the effective Bohr radius, ao = eh 2 / '/ie 2 , where /x 
is the electronic effective mass, and e is the dielectric 
constant of the semiconductor medium. For instance, 
in the particular case of GaAs-based semiconductors, 
fi = 0.067m e , and e = 13.18 . Thus, for a GaAs host, 
the effective Rydberg is numerically Ry = 5.2 meV, and 
the effective Bohr radius ao = 104 A. 

The rest of this paper is organized as follows. In 
section 2, the theory of electron energy states spectrum 
and wave functions is developed for the MLQAD with 
hydrogenic donor impurity at the center. In theoretical 
calculations, effective mass approximation and spher- 
ically symmetric are employed. In this sense, the ex- 
act solution of the Schrodingcr equation is expressed in 
terms of the Whittaker functions. In section 3 the en- 
ergy spectrum, radial probability distribution and os- 
cillator strength of the intersubband transitions are ob- 
tained for different cases. Finally, our conclusions are 
presented in the last section. 

2 Theory and formulation 

In the effective mass approximation, the Hamiltonian 
for an on-center hydrogenic impurity in MLQADs is 
given by 

H=-^-V 2 + U(r) (1) 

where fi is the electronic effective mass in the semicon- 
ductor medium and U(r) represents the total potential 
(coulomb and CP) as following 

f Vi - r <a 

U(r) = lv 2 - |, a < r < b (2) 

r>b 

where a is the core radius and b is the total antidot 
(core plus shell) radius, therefore b — a is the thickness 
of the shell. V\ and V% are the positive CP inside the 
core and shell media respectively. The relation between 
the value of V\ and V2 is V\ > V2, this condition ensures 
that the electron reminds bind in several lowest state 
energies at least. Furthermore, the value of V3 is always 
zero (i.e., the bulk material have no CP) but writing 
V3 helps us to be able for writing the next equations in 
the compact form as you will see. 

In the three dimensions spherical coordinate, the 
Schrodinger eigenvalue equation is Hip(r, 9, <ft) — Eip(r, 9, 
where the eigenfunctions are of the separable form ip{r, 9, 
R(r)Y (9 , 4>) . For the r-dependent (spherically symmet- 
ric) potentials, the angular dependence of the wave func- 
tion Y(9, </>), is given by the spherical harmonics [25] . 
In particular, the angular part is unaffected by the core 
and total radius, thus we ignore this part. The radial 
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part is affected by the core and totai radius and the 
antidot CPs, so in the following it is considered. The 
radial part of the Schrodinger equation in the spherical 
coordinate for three regions are given by the following 
compact relation 



electron energies and normalized wave functions the os- 
cillator strength of the intersubband transition between 
the states ipnim and ipn'i'm' in the dipole approxima- 
tion, with the radiation of linearly polarized light along 
z-axis, are given by the following expression [7 



2/i dr 2 r dr r 2 

where i = 1, 2 and 3 are corresponding to r < a, a < 
r < b and r > b respectively. Using the convenient 
parameters 

2/ie 2 



0(3) 



2 _ -8fi(E - Vj) 



eh A dj 



(4) 



and further writing W(£i) = £iRy'(r), then the above 
Schrodinger equation conveniently converts to the fol- 
lowing whittaker equation 



The general solution is the whittaker function 



(6) 



where Cu and C21 are constants to be specified by con- 
tinuity, asymptotic and normalization conditions. The 
asymptotic behavior dictates that the wave function 
must be finite and zero in origin neighborhood and far 
away from origin, respectively. These conditions lead to 
C21 = C13 = (since W Ai)J+ i(&) and M A<) , + i(&) di- 
verge at the origin and infinity, respectively). Now we 
can write the wave functions in three regions as 



= C n 
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r < a 



(7) 



, 2) c 13 M AaiJ+i (6)+c 22 w AaiJ+i (e 2 ) 

B> > = a < r < b 
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> b 



(8) 
(9) 



The boundary conditions of the radial wave functions 
and their densities of current continuity at the inter- 
faces r = a and r = b require 

dr ctr 

ctr dr 
By using the above two equations, one can obtain the 
eigenvalue E of electron in MLQAD with donor impu- 
rity in the center. 

The oscillator strength is a very important physi- 
cal quantity in the study of optical properties which 
related to the electronic transition. Using the obtained 
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n>r 3 dr 



Y lm Y; m , cosOdQ 



(12) 



where, AEfi denote the energy difference between the 
final and the initial state energy and |(/|.z|i)| 2 is the 
transition matrix element that in second line it has 
been expressed in terms of normalized wave functions 
and spherical harmonics. The selection rules for such 
transitions are determined by integrating over angular 
variables. With the substitution cos 9 = (V47r/3)Yi )0 , 
it can be readily shown that the probability of tran- 
sition is unequal to zero only when \l — l'\ — 1 and 
m = m' , where I is the orbital, and m is the magnetic 
quantum number. By using of this information and the 
layered nature of MLQAD, one can rewrite the oscilla- 
tor strength as 



fn 



lm—n' V m 



, = 2 4ae 



fi 



Y lm Y{; m , cos 9df2 



R^R^dr- 



p(3) d(3) 



3 Results and discussion 

3.1 Energy spectrum 

Here, the binding energies of the ground state and the 
first two excited states of a hydrogenic impurity located 
at the centre of the MLQAD are calculated as a function 
of core radius for different shell thicknesses and CPs. To 
make a comparison, we focus on the special case of the 
same CPs (i.e., V\ = V2 = V). In fig.ffl the energy spec- 
trum for Vi = V2 = 2Ry as a function of antidot radius 
are plotted. As we expected, the result are identical to 
the SLQAD that has been studied in [TIE]. It is obvi- 
ous from fig. [1] that the orbital energy degeneracy in 2s 
and 2p states, is removed. The energies for a free hydro- 
gen atom is — l/n 2 Ry, where n is the principle quan- 
tum number, so the ground state (Is) and the first two 
excited state (2s and 2p) energies are equal to —lRy 
and —§.2bRy respectively. From fig-U one can see that 
the energy curves will approach to the energy values of 
a free hydrogen atom in each state when the antidot 
radius becomes smaller and smaller. Although, when 
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Fig. 1 The binding energies of the ground, first and sec- 
ond excited state of the MLQAD hydrogenic impurity for the 
same confining core and shell potential (V± = Vb = 2Ry) as a 
function of antidot radius for (a) finite CPs (b) infinite CPs. 

the radius of impurity increases, the confinement effect 
enhances the state energies and at very large radii, a 
bounded electron to an impurity never sees the surface 
of the antidot potential and behaves as a free electron. 
In the case of a MLQAD with different CPs {V x ^V 2 ), 
the ground state energy is plotted for two different cases 
in fig. [2] In fig. [2] (a) the ground state energy vs total 
antidot radius is plotted when the shell CP and the 
core radius are fixed as V% = 2Ry and a = ao whiles a 
core CP V\ has three distinct values 5, 10 and ooRy. 
In fig. [5] (b) the variation of the ground state energy 
with respect to the core radius is investigated for fixed 
values of CPs and shell thickness. From fig. [5] (a), it 
is clear that the binding energies are not influenced 
by the core CP value when the total radius is larger 
than about 2.5ao but they are dependent to the core 
CP value when the total radius is less than 2.5ao- More 
precise, in the range of ao < b < 2.5a , a higher core CP 
leads to a larger binding energy, fig. [5] (b) shows that in 
the small core radii, the smaller shell thickness leads to 
the more bounded electron, while in large enough core 
radii (about a > 8ao) the results are independent on the 
shell thickness and the electron is no strongly bind in 
this case. Also, when the shell thickness becomes very 
small the ground state energy curve approaches to the 
SLQAD case (see fig. Q] solid curve). 

3.2 Radial probability distribution 

In figs. [3] (a) and (b) the radial probability distribution 
of the electron in its ground state, Pi s , for a couple of 
core and shell CPs and for three different total antidot 
radii, are plotted. The results for the radial probabil- 



ity distribution, clearly indicate that by increasing the 
total radius, the radial probability function becomes 
broader and expands to longer distances from the ori- 
gin, in this sense, the electron becomes more delocal- 
ized. As the radial probability distribution gets broader, 
its height decreases, in accordance with the conserva- 
tion of probability. By comparing the mutual curves 
(solid by solid and so on) we found that the core CP 
values have influenced mainly on the radial probability 
distribution of antidot with the small total radius (i.e., 
b = 1.2ao or solid curves) and two other curves (dashed 
and dotted curve) are closely similar for both core CPs 
V\ = 5Ry and V\ = oo. The radial probability distribu- 
tion of the ground and second excited state for various 
shell thicknesses are computed in figs. 0] (a) and (b) 
respectively. These graphes indicate that by increasing 
the shell thickness the radial probability density func- 
tion expanded to farther places from the origin. In the 
limiting case (b — a) — >• the result approaches to the 
SLQAD case as expected. 



3.3 Oscillator strength 

The oscillator strength of an electron bound to a donor 
impurity for transition from the ground state to the 
first excited state, and for various shell thicknesses with 
finite and infinite core CP, are plotted as a function of 
core radius in figs. [5] (a) and (b). As seen from figs, 
the oscillator strength curves undergo a peak in given 
core radii due to the big energy difference and good 
overlapping of the corresponding wave functions. When 
the value of shell thickness increases, the place of peak 
shifts toward smaller core radius and when the shell 
thickness becomes small (b — a — >• 0), the result (solid 
curve) approaches to the SLQAD case. By comparing 
the figs. [5] (a) and (b) one can find that by increasing the 
core potential (from 5Ry to oo) the difference among 
different shell thicknesses becomes smaller. 

In fig. |6l the variation of oscillator strength with re- 
spect to the total antidot radius for several core CPs 
have given when the core radius and the shell CP are 
fixed as a = 0.5ao and Vi = 2Ry respectively. Accord- 
ing to the normalization condition, the summation of 
all allowed transition probabilities is 1. From fig. H](a), 
one can easily see that at the given core radius (about 
a = 2ao), /2 io- ioo displays a maximum and obtains 
the major portion of I (larger than 0.95), so in this 
given radius the other transition probabilities have a 
very small portion of I and tend toward 0. Fig. [6] (b) 
shows that at about a = 2a>o that /210-100 displays a 
maximum, /310-100 displays a minimum and tends to- 
ward 0. However, for large enough total antidot radii, 
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Fig. 3 The radial probability distribution of the ground state (-Pio) for several total radii with a = Oq by choosing (a) 
Vi = 5Ry and V 2 = 2Ry (b) Vj_ = oo and V 2 = 2Ry. 




Fig. 4 The radial probability distribution for various shell thicknesses b — a — 0.2,0.5, lao with a — ao, V± = 5Ry and 
V2 = 2Ry for (a) the ground state, Pio (b) the second excited state, P20- 
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Total radius (a ) Total radius (a ) 

Fig. 6 Oscillator strength from the ground to (a) the first allowed excited state (/210-100) (t>) the second allowed excited 
state (/310-100) f° r several confining core potentials V\ = 5, 10, ooRy with constant shell potential V2 = 2Ry and a = 0.5ao 
as a function of the total antidot radius. 



the two transition probabilities tend to become compa- 
rable. 



4 conclusion 

The exact solution of Schrodinger equation within the 
effective mass approximation are obtain for the electron 
in the MLQAD. It is found that only in the small total 
antidot radius a higher core CP leads to a larger binding 
energy. Also, it is shown that in the small core radii, the 
smaller shell thickness leads to the more bounded elec- 
tron, while in large enough core radii (about a > 8ao) 
the results are independent on the shell thickness. The 
oscillator strength for the intersubband electric transi- 
tion from the ground state to the two first allowed ex- 



cited state calculated by using the obtained spectrum 
and the corresponding wave functions. It found that for 
a fixed value of the core radius and shell potential by 
varying the core CP value, /210-100 displays a maxi- 
mum larger than 0.95 at the given core radius (about 
a = 2ao) so the transitions involving the second allowed 
state, /310-ioOj can be neglected. But, in particular, for 
large enough total antidot radii, /210-100 and /310-100 
become comparable. This indicates that the transitions 
involving the second allowed state, are no longer neg- 
ligible when investigating the optical properties of the 
hydrogenic MLQAD with large enough total antidot 
radii, of the order of several hundred nanometers. 
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